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Predicting the Influence of Plate Geometry on the Eddy Current Pendulum 


Catherine Weigel, Jeremy Wachter, Paul Wagoner, and Timothy J. Atherton* 
Department of Physics and Astronomy, Tufts University, Medford, MA 02155, USA 

We quantitatively analyze a familiar classroom demonstration. Van Waltenhofen’s eddy current 
pendulum, to predict the damping effect for a variety of plate geometries from first principles. 
Results from conformal mapping, finite element simulations and a simplified model suitable for 
introductory classes are compared with experiments. 


I. INTRODUCTION 

Eddy currents are induced electric currents in a con¬ 
ducting material that result when either the object moves 
through a nonuniform magnetic field or is stationary but 
subject to a time-changing magnetic field. They were 
first observed by Frangois Arago[l] but it was not un¬ 
til Michael Faraday’s discovery of induction[2] that the 
mechanism of this phenomenon was understood. A pi¬ 
oneering experiment to quantify the heat dissipated by 
the eddy currents, and hence connect electromagnetism 
to thermodynamics, is von Waltenhofen’s pendulum[3] 
[figure 1(a)— (c)]. The bob of the pendulum is a non fer¬ 
romagnetic conducting plate which swings between the 
poles of an electromagnet, generating circulating eddy 
currents. The dissipative eddy currents tend to oppose 
the motion of the bob, and hence damp the pendulum. 
The experiment is an example of magnetic braking: when 
the magnet is turned off, the pendulum oscillates as a 
classic pendulum, damped only by air resistance and 
frictional forces in the hinge; when the magnetic field 
is switched on, the pendulum comes to a rapid stop. 

Von Waltenhofen’s pendulum is a popular demonstra¬ 
tion in high school[4] and introductory undergraduate [5] 
classes because it provides a striking visual illustration of 
the invisible eddy currents. It is easy to set up, and gives 
the instructor a chance to discuss real-world applications 
such as magnetic braking in roller coasters, hybrid cars 
and the detection of flaws in conductive materials. More¬ 
over, by using different shaped plates as bobs, some of 
the characteristics of eddy currents can be explored: if 
the plate is now substituted with one where the slits do 
not reach the edges [figure 1(d) (ii)[, the damping is only 
mildly altered. If a solid plate is replaced by a comb-like 
plate with slits cut into it [figure l(d)(iii)-(iv)[ the damp¬ 
ing effect is greatly attenuated. Students infer from this 
that the eddy currents must be localized to some region 
of the plate. 

Typically, however, the relationship between the shape 
of the plate and the damping is considered only qualita¬ 
tively due to the complexity of the calculations required, 
which require advanced mathematical techniques intro¬ 
duced in higher-level electromagnetism classes. In this 
paper, therefore, we apply some of these strategies to 
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predict the damping behavior of conducting plates with 
different geometries. We investigate four plate shapes, 
illustrated in figure 1(d): (i) a square plate, (ii) a plate 
with holes that do not reach the edge, (iii) a plate with 
with two slots, (iv) another with four slots. The present 
work was primarily performed during a project-based 
graduate electrodynamics class with the objective of giv¬ 
ing students a real, tractable example of the utility of 
their calculations. It was also an opportunity to address 
the challenges of connecting theory to experiment, even 
where the fundamental theory is uncontroversial. 

The paper is organized as follows: in section II, the 
form of the damping for the present experiment is de¬ 
rived first in terms of a simplified model suitable for 
undergraduate classes and the effect of these damping 
terms on the simple pendulum is explored. In section 
III, we present a more sophisticated model derived from 
Maxwell’s equations in the magneto-quasistatic approx¬ 
imation. Within this framework, the rate of power dis¬ 
sipation for each plate is estimated using a conformal 
mapping technique and a finite element solution. In sec¬ 
tion IV, the motion of the pendulum with each of the 
four plates is determined experimentally and compared 
with the predictions of each of the theoretical models in 
section V. In section VI, the pedagogical context of this 
work is discussed. Brief conclusions are presented in sec¬ 
tion VII. 


II. SIMPLIFIED MODEL 

The damped simple pendulum is, in the presence of 
dissipative forces, described by the equation, 

■■ 2 dD ^ 

q + io^qP —Q ( 1 ) 

where D{q,q) is the Rayleigh dissipation function, wq 
is the natural frequency of the pendulum and q{t) is a 
generalized coordinate. In the absence of the magnetic 
field, a Stokes-like drag is assumed, 

D stokes = [iif, (2) 

arising from friction in the bearing and air resistance. In 
this section and the next, the form of the dissipation func¬ 
tion in the presence of the magnetic field is derived, first 
for a simplified model suitable for introductory classes 
and second a more physically correct model based on 
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Figure 1: Schematic of Von Waltenhofen’s pendulum, (a) 
An aluminium plate swings between the pole pieces of an 
electromagnet. As it passes through the fringing field, eddy 
currents are induced in the plate that dissipate energy and 
damp the pendulum, (b) As the plate passes through the 
center of the pole pieces, the rate of change in flux through the 
plate vanishes and so the eddy current momentarily vanishes, 
(c) As the pendulum swings back, the eddy current reverses 
direction, (d) Various plate geometries under investigation in 
the present work (to scale). 


the magneto-quasistatic approximation. The remaining 
theory section considers the effect of these terms on the 
solutions of (1). 

To develop the simplified model, we make the gross 
approximation that the eddy currents are completely lo¬ 
calized to the boundary of the plate. Then, the plate 
may be replaced by a thin closed loop of wire of equiva¬ 
lent shape with area A, length L, resistivity p and with a 
normal unit vector n. The magnetic flux passing through 
the loop is, 


$ = 


B • dA = (B) A 


( 3 ) 


where (B) is the mean value of the magnetic field over the 
plate when it is located at a particular position q. The 
emf induced around the loop is, by the Faraday-Lenz law 
and using the chain rule. 
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9$ 

W 


. d (B • n) 


4 - 


( 4 ) 


The power dissipated in the loop is precisely the dissipa¬ 
tion function desired and readily calculated, 

D = V^R, (5) 


which, using R = pL, becomes, 


/A2\ 1 /d(B-n) 
V L j p V dq 


2 


( 6 ) 


where the highlighted term represents the “geometric fac¬ 
tor” arising from the shape of the plate. The dissipation 
function (6) has the form of a position-dependent Stokes 
drag. 

It is now necessary to solve the pendulum equation 
including the new dissipation terms (6), 


q -f- uJq q -f- 2 


/d(B-n) 
Lp \ dq 



( 7 ) 


for which the form of the magnetic field must be known. 
Here, the coordinate system is oriented as shown in figure 
1 such that the plate normal n is parallel to the z axis, 
so only is required; the origin is chosen to lie in the 
middle of the pole pieces. An approximation made in 
other studies[4, 5] is that the magnetic field is constant 
between the pole pieces and zero everywhere else, i.e. 
that, 




Bq Bpp / 2 ^ X < Lpp! 2 

0 X < —Lpp/2, X > Lpp/2 


( 8 ) 


where Lpp is the width of the pole pieces. If the plate is 
square of side L, the average value of the field over the 
surface of the plate when it is centered at q along the 
x-axis may be determined, 


l‘2 


(B • n) = L 


/ q-L/2 


Bz{x') dx'. 


( 9 ) 


For the ansatz (8), the flux through the plate increases 
(decreases) linearly as it enters (leaves) the field. De¬ 
pending on the size of the plate, and the amplitude of 
the motion, the plate may also experience portions of 
the motion where there is no change in flux. If the 
plate is the same size as the pole pieces L = Lpp, and 
never fully leaves the field, i.e. |x| < {Lpp -|- T)/2 , then 
(B • n) ^ B^q and the classical damped pendulum is re¬ 
covered, but with a modified damping constant. 


X -t- COnX 2 


42 n2 


Lp 


X = 0. 


( 10 ) 


In order to gauge the validity of the (8) for the pendu¬ 
lum we calculated the spatial dependence of the magnetic 
field due to the pole pieces in the magnetostatic formula¬ 
tion, using a scalar potential function defined B = — 
and subject to Laplace’s equation = 0. This was 
solved using the finite element software FlexPDE for 
square pole pieces of width 11cm and separation 3.5cm. 
Dirichlet boundary conditions were imposed on the pole 
pieces. A plot of B^ along the horizontal centerline of 
the midplane of the pole pieces is displayed in fig. 2(a), 
showing as expected that the field is constant between 
the pole pieces with the fringing field decaying rapidly 
outside. From this calculated profile, the flux J B • dA 
through a square plate placed at different horizontal po¬ 
sitions was also calculated; the result is displayed in fig. 
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Figure 2: (a) Field strength as a function of distance x 

through the center of the pole pieces, (b) Magnetic flux 
through a square plate as a function of the position of the 
plate X. 


2(b). Notice that the flux depends linearly on position 
for 1 ^ |a;| < lOcm, suggesting that the equation of mo¬ 
tion (10) is likely to be valid for moderate amplitudes; 
however around a; = 0 the flux depends quadratically on 
X, suggesting that (10) is likely to be inaccurate for small 
amplitudes. 


III. MAGNETO-QUASISTATIC MODEL 

Since the eddy current is not strictly confined to the 
edges of the plate, it is necessary to develop a more so¬ 
phisticated model by solving Maxwell’s equations 

niD 

V-E = 0, V-B = 0, VxH = J, VxE-h—=0 (11) 

ot 

together with the continuity equation 

V.J + |=0 (12) 

to determine the true current density in the plate. 
Here, the magneto-quasistatic approximation [6] has been 
made, i.e. the Maxwell correction to Ampere’s law is ne¬ 
glected. This is justified because the timescale of varia¬ 
tion of the magnetic field over the surface of the plate is 
sufficiently long that the eddy current produced does not 
significantly alter the magnetic field in turn. 

To perform the analysis, consider a thin metal plate 
of uniform conductivity a and thickness d and where the 
shape of the plate is arbitrary. It is assumed that Ohm’s 
law holds for the plate, so 

J = crE. (13) 

In the frame of the laboratory, since the pendulum is 
rigid, the position of the plate may be described by the 


single generalized coordinate g(t) in eq. 1. As the pendu¬ 
lum swings, the plate passes through regions of different 
magnetic field strength. Viewed from the point of view 
of the plate, the magnetic field across its surface changes 
as a function of time. This may be described by defin¬ 
ing a local coordinate system {x',y',z') relative to the 
plate such that the plate lies in the x' — y' plane, and 
considering the magnetic field in this frame, 

B = B(cr', 2 /',g(t)). (14) 

We note that the reference frame {x', y', z') is not an in¬ 
ertial frame as the pendulum undergoes periodic acceler¬ 
ation during the course of its motion; however we neglect 
these effects as the pendulum velocity rq <C c. By Fara¬ 
day’s law, an electric field and hence an eddy current is 
induced in the plate: if d <C 5 the skin depth of the ma¬ 
terial it may be assumed that E is oriented in the x — y 
plane and does not vary significantly over the thickness 
of the plate; hence the problem is quasi-two-dimensional. 

From equations (11) and (13), it is a standard textbook 
problem to show that the quantities E, B, A and J all 
obey the diffusion equation, 

WJ=4. (15) 

Since p = 0, the continuity equation implies that J must 
have no divergence and hence can be written using a 
stream function, 

d(j){x',y') d(j){x',y’) 

■'» = — 

By inserting (16) into (15), together with the given ap¬ 
plied Magnetic held into Maxwell’s equations, it may be 
shown that the stream function obeys Poisson’s equation 

V^(/)(a;', y') = ^ q = BoaqC{x', q) (17) 

where V here is the 2D Laplacian, Bq is the magnetic held 
strength and we defined a source term C(x, q) = /Bq. 
Since J must lie tangentially to the edge of the plate, 4> 
obeys a Dirichlet boundary condition (/> = 0 on all edges 
of the plate. From the solution for (p, the instantaneous 
rate of power dissipation by the plate is then calculated, 

D=^l\3fdA, (18) 

where the integral is to be taken over the plate. Inserting 
the stream function (16) into (18), integrating by parts 
and using (17), we obtain 


J vp-vp dA, 

(19) 


(20) 

dqBo [ (AC(x',g) dA. 

(21) 
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A formal solution to (17) may be constructed using the 
Green function for the Laplacian, 


Inserting this into (21) and rearranging we obtain the 
dissipation function, 


= aqBoGC = <jqBo j G(x',x")C(x",g)dA". (22) 


D{q,q) = -dBoq j aB^q ^ G(x', x")C(x", g) dA" 


C(x') dA' 


-daB^q^ j J ^(x', g)<^(x", g)G(x', x") dA'dA". 


(23) 

(24) 


Since the source function C(x, q) depends on the position 
of the pendulum q, this has the form 

D{q,q)^BlP{q)e (25) 

where we define a local dissipation rate P{q). To de¬ 
termine an effective damping coefficient, the calculations 
proceed as follows: in subsection III A the form of ()(x, q) 
is determined, then the positional dependence of the drag 
term is evaluated using a conformal mapping technique 
in subsection IIIB, and also by solving (17) with FlexPDE 
in subsection IIIC; finally, the effect of the position de¬ 
pendent damping on the motion of the pendulum is con¬ 
sidered in the subsection HID. 


A. Evaluating the Fringe field 

Before evaluating the power dissipation, it is first nec¬ 
essary to determine the functional form of the magnetic 
fringe field that the plate experiences as it moves between 
the pole pieces and hence the function C(3:)- This could 
also be obtained from a finite element calculation as was 
done for the simplified model in section (II) . We made the 
approximation that only the fringe fields from the sides 
of the pole pieces are significant in the damping, and 
the vertical extent of the pole pieces are neglected; the 
formulation above remains sufficiently general that these 
details could easily be included in a more sophisticated 
calculation of C{x). 

Within this approximation, the problem is two dimen¬ 
sional and hence can be solved using conformal mapping, 
an important technique for solving electrostatics prob¬ 
lems. For a full discussion we refer the reader to [7] and 
[8]. Briefly, however: A map from a domain to its image 
is called conformal if angles measured locally around a 
point in the source domain are preserved under the map¬ 
ping. A key property of such maps is that if a function 
that is a solution to Laplace’s equation is constructed on 
some domain, then the projection of that function under 
the action of the mapping will also obey Laplace’s equa¬ 
tion. Boundary value electrostatics problems in two di¬ 
mensions can hence often be solved using this method by 



Re(w) Re(z) 


Figure 3: Conformal mapping z{w) from (a) a portion of the 
strip 0 < A(w) < 1 to (b) the depicted region in the 2 plane. 
Places where the integrand in (26) vanishes are indicated by 
red circles; these introduce corners into the map. Contours of 
constant K(w) (dashed lines) and Q{w) (solid lines) and their 
projection are indicated corresponding to equipotentials and 
field lines respectively. Inset: Plots of the field profile B^ix) 
and its spatial derivative (x) calculated from the conformal 
mapping. 

constructing a solution to Laplace’s equation with appro¬ 
priate boundary conditions in some simple domain and 
then finding a conformal map to the domain of inter¬ 
est. The powerful Riemann mapping theorem guarantees 
such a map exists if the domains are simply connected; 
mappings between multiply-connected domains may also 
exist if certain compatibility criteria are met [8]. A rich 
source of such maps are analytic complex functions z(w) 
that map a domain in the complex plane w = u + iv to 
an image in the plane z = x(u, v)+iy(u, v). The required 
mapping to find the fringe field is, 

nW 

z{w) = / exp(—7r^/2) [sinh (tt^)]^'^^ d^. (26) 

Jo 

As shown in figure 3, this maps the strip defined by 0 < 
3(?u) < 1 in the w = u+iv plane onto the depicted region 
in the z = x + iy plane; the resulting shape resembles an 
overhead view of the edge of two parallel pole pieces. The 
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integrand in (26) vanishes at ^ = 0 and ^ = i- The effect 
of these poles is to introduce corners in the map; the 
angle of these is controlled by the strength of the pole, 
here 1/2. The factor of exp(—7r^/2) is required to ensure 
that the boundary of the 3?(w) < 0 part of the strip 
remains horizontal while the ^(w) > 0 half is mapped to 
the vertical boundaries. The integral can be performed 
analytically, 


z(w) 




e'^“' arcsin (e^*") \ 

Vl - J 


V^sinh/TTw), 


(27) 

where the overall prefactor is chosen to ensure that the 
3?(?u) < 0 half of the strip is mapped onto 0 < 3(z) < 1. 

The magnetic field strength can be determined using 
this map as follows: in the space between the pole pieces 
V-B = 0 and hence the magnetic field can be constructed 
from a scalar potential B = —Vip subject to Laplace’s 
equation V^ip = 0. On the boundary, B must be per¬ 
pendicular to the pole pieces so these must be surfaces of 
constant ip. In the w-plane, the system resembles a paral¬ 
lel plate capacitor with the solution ip = Av -I- B where A 
and B are arbitrary constants. Since Laplace’s equation 
is invariant under the conformal mapping, this solution 
remains valid when projected into the z-plane: the lines 
of constant v are equipotentials and lines of constant u 
are the field lines as shown in figure 3. 

The magnetic field can be found from B = 
—A {dxV, dyv) ; this can be expressed in complex form 
in the w-plane as B{w) = —A/z’{w)\l\. Evaluating z{w) 
and B{w) along the center of the strip w = t + i/2 yields 
By along the midplane; the x derivative of By required 
for the pendulum motion is evaluated from these using 
the chain rule. 


X 


V 

dx 


- [e’^Vl + e-2’^* 

TT L 

—-7rsech(7rt)^. 


2arccosh(e’"^*) 


(28) 


Having tabulated values of B'y{x) using the above im¬ 
plicit expressions, a final form for (p{x) was assembled 
from two appropriately translated and rescaled copies of 
By{x) as shown in the inset of figure 3; this was used in 
subsequent calculations as an InterpolatingFunction 
in Mathematica, or exported as a table for use in 
FlexPDE. For the plates and pole pieces considered, it 
is actually the case that the plates were slightly larger 
than the width of the pole pieces; while this could easily 
be accommodated by a simple rescaling of the solution 
(28), it was found in practice not to affect the dissipation 
enough to affect the quality of the fit. 


B. Conformal Mapping Solution 

Conformal mapping can be also used to evaluate the 
power dissipation of the plate as it passes through the 





Figure 4: (a) Map from the unit circle to a paddle with a single 
slit of variable geometry, (b) Map from a rectangular domain 
to the slit. The position of poles in the canonical domains is 
indicated with red circles; these become vertices in the image 
domain, (c) Two maps relevant to the experiment: i) map 
from a circle to a square; ii) map from a rectangular domain 
to a square with two slits—note significant crowding in the 
central spar. 


fringing field. The strategy here is to find a mapping from 
a simple canonical domain on which the Green function is 
known to one that approximates the shape of the plates; 
then the integral (24) can be computed. As will be seen, 
the choice of canonical domain is important to render the 
problem tractable. We therefore examined two choices: 
first, the mapping from the unit circle (w-plane) to a 
polygon of interest (z-plane) defined by an ordered set of 
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vertices is given by the Schwarz-Christoffel formula[8], 

/ W _ 

( 29 ) 


where at is the exterior turning angle of the polygon at 
vertex i and the are the points on the unit circle 
that are mapped to the vertices of the polygon The 
constants A and C are translation and scaling factors re¬ 
spectively. While the set of parameters {ai} can be de¬ 
termined directly from the polygon, the positions of the 
parameters {w®} must be determined numerically[8]. A 
and C are then determined to rotate and scale the result 
onto the polygon of interest. Some representative results 
for a single slit are shown in fig. 4(a). Due to the sym¬ 
metry of the problem, the position of three poles must 
be determined numerically with the remaining vertices 
provided by reflection. 

The second choice of canonical domain is a rectangular 
portion of the semi-infinite strip, —W < 3fi(?u) < W , 
0 < 3(w) < 1. A mapping onto the shape of interest is 
given by. 


z = A + C 



'TT 

.2 





(30) 


where similarly are the position of the poles which 
produce the required corners in the image domain and 
at are the turning angles; this formula is a more general 
version of equation (26) used to determine the fringing 
field. Results for the single slit are shown in fig. 4(b). 
Note that this map is approximate only: the left and 
right sides of the rectangle are mapped to the bottom of 
the paddle—the side that the slit is on—but their image 
under (30) is not a straight line. Hence, only two poles 
{w®}—mapped to the interior vertices of the slit—need 
be determined numerically; the third parameter is the 
width of the rectangle W. Despite the approximate na¬ 
ture of the map, the distortion of the lower boundary is 
very small. These parameters were found by minimizing 
an objective function that was constructed from the L 2 
norm of the difference between the ratios of the sides of 
the polygon and their desired values; a notebook to do so 
is provided as Supplementary Material. Robust numer¬ 
ical codes including SCPACK[9] and the SC Toolbox[10] 
are available for more complex shapes. 

Comparing the two choices of canonical domain in fig¬ 
ure 4, an important difference is apparent. The poles in 
the maps from the unit circle tend to lie very close to 
one another, becoming in some cases indistinguishable in 
the figure, though they remain numerically distinct. As a 
consequence of this phenomenon, known as crowding[S\, 
very small regions adjacent to these crowded poles are 
mapped to large regions in the image domain. Crowding 
is problematic, because it makes numerically finding the 
challenging; it will also cause problems in evaluating 


the power dissipation. Notice that using the rectangular 
source domain eliminates the crowding problem as the 
poles remain well-separated. 

In fig. 4(c), we show two maps relevant to the exper¬ 
iment. First, the map from the circle to the square is 
famously 


An example with two slits that closely matches the ex¬ 
perimental plate is shown in fig. 4(c) (ii), constructed by 
introducing additional poles into one side of the one-slit 
rectangle maps above. Unfortunately, the crowding phe¬ 
nomenon is visible in this map: a very small piece of the 
rectangle is mapped to the central spar of the plate. As 
will be seen later, this will reduce the accuracy of the 
power dissipation calculation when this spar coincides 
with the edge of the pole pieces. While it is possible to 
find maps with crowding localized to different portions of 
the plate[8], it cannot be alleviated entirely if more than 
one slit is present. Hence, we did not pursue maps for 
plates with more than 2 slits. 

Turning, then, to the issue of evaluating the power 
dissipated by the induced current, we initially tried to 
use the eigenfunction expansion of the Green function, 

V'n(x)V'n(x') 

,X j — > ---, 

which, inserted into the integral in (24) allows 
obtain, 

J j C(x',g)C(x",g)G(x',x") dA' dA" 

[/C(x,g)V’n(x) dA]^ 

n 

The 2D integral in (33) can then be performed in the 
w-plane using the normalized eigenfunctions of the unit 
disk and the mapping z(w), 


(32) 

one to 

(33) 

(34) 


J C(x)^/’n(x) dA = y Ciziw),q)ipniw)\z'{w)\^ dw. 

(35) 

Unfortunately, we found that the series (34) converges 
very slowly except for the square map (31). The problem 
is further compounded for the unit circle by the fact that 
the integrand in (35) is very sharply peaked and localized 
to the boundary in many cases due to crowding. For these 
reasons, we abandoned this approach and the unit circle 
maps as viable means to evaluate the dissipation. 

The power dissipation integral (24) can be re-expressed 
in the w-plane as. 
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D (X J J (!^{z{w'),q)(^{z(w''),q)G{w\w'')\z'{w')f \z'{w'')f dw' dw". (36) 


In order to evaluate this integral, one requires a numeri¬ 
cally efficient representation of the Green function. Un¬ 
fortunately, the Green’s function for a rectangle cannot 
be expressed in closed form, requiring expensive evalua¬ 
tion of elliptic integrals. Recently, however, renewed in¬ 
terest in Green function methods has yielded rapidly con¬ 
verging series representations[ll]. Briefly, the strategy 
is to start from the eigenfunction expansion of G(w,w') 
(32) on the rectangle and separate the double summation 


into two pieces: one that can be summed analytically and 
a second that remains to be performed numerically. For¬ 
tuitously, the analytical part contains the inherent sin¬ 
gularity; the remaining numerical sum is a smooth func¬ 
tion requiring few terms to converge to machine precision. 
Equation (7) of [11] presents an expression for the green 
function in a square; we derived an equivalent expression 
for the rectangle —W < ?ft(w) < W , 0 < Q(w) < 1 of 
interest. 


G{w,w') = 


^ In. 

1 - 

p _ qW+w' 

o 

, V 

sinh (riTTu) sinh (jittu') 

27r 

1 - 

2 _ 

TT 

n—1 

j^g 2 n-nW ( 2 n 7 rkF) 


sin (riTTv) sin (mrv') 


(37) 


where w = u + iv and w' = u' + iv' consistent with our 
earlier definitions. 

Inserting (37), (30) and ({x) derived in the previous 
section into (36), we evaluated the local power dissipa¬ 
tion rate P{q) as a function of q for a variety of plate 
geometries. A notebook to perform these calculations 
is presented as Supplementary Material; we defer con¬ 
sideration of the results until the subsequent section, in 
which the same function is calculated by finite element 
simulations. 


C. Finite Element Simulations 

Finite element analysis is a commonly used numerical 
technique for solving partial differential equations in sys¬ 
tems too complicated for analytical solution. Here, the 
PDE of interest is, as shown in the previous sections, 
Poisson’s equation with a spatially dependent source 
term (17). This must be solved to obtain the stream 
function , y') for each of the paddle shapes of interest 
and at each point q in the motion of the paddle. Suit¬ 
able input files for FlexPDE for each plate are included as 
Supplementary Material. Each script contains field vari¬ 
able and parameter definitions, a definition of the PDE, 
a description of the appropriate domain boundary and 
specification of the boundary conditions, and a list of 
the outputs desired, i.e. (j), J and the value of the in¬ 
tegral (24). A relative error of < 10“^ in the absolute 
(p was requested and used by FlexPDE to guide adaptive 
refinement. 

For each plate, the boundary condition on the exte¬ 
rior is simply (p = 0, which forces the stream function 


to lie tangentially to the edge of the plate. The plate 
with holes requires special treatment: around each of 
the holes p should be constant, i.e. the hole should be 
an equipotential surface, but the correct value is not a 
priori known. A way to obtain the correct solution is 
to exploit an electrostatic analogy, modifying Poisson’s 
equation by introducing a permittivity e, 

V ■ieVP) = ax',q). (38) 

The computational domain is extended so that p is also 
defined inside the holes. Different values of e are used for 
the interior of the plate Cpiate and the holes tholes- As the 
ratio tholes/tpiate oo, the solution converges to the de¬ 
sired solution where p constant inside the holes. For the 
purposes of calculation, tpiate is chosen arbitrarily to be 1 
and tholes is successively increased; the value of (24) con¬ 
verged to a relative error of 10“^ for tholes > 10^. Plots 
of the stream function p and the corresponding current 
densities J resulting from the calculations are displayed 
in figure 5, together with a plot of the dissipation rate 

P{q)- 

It is now possible to compare the Gonformal Map¬ 
ping and Finite Element methods for this problem. To 
do so, we used both methods to calculate the dissi¬ 
pation rate P{q) for the plate geometries displayed in 
fig. 4(b) and (c), i.e. those for which we were able 
to And appropropriate mappings. To perform the inte¬ 
gral (36), we used Mathematica’s NIntegrate with the 
AdaptiveMonteCarlo method and 10® evaluation points. 
The calculation for each plate took ^ 5 minutes for 
FlexPDE and ^ 20 minutes for the conformal mapping 
approach on the same computer (Apple Macbook Air 11” 
Late 2014). Results for the various plates are shown in 













Figure 5: (a) Dissipation rate as a function of position and 
(b) i-iv corresponding current density profiles obtained nu¬ 
merically for the various shapes from FlexPDE. 


fig. (6). The agreement is very good for the single slit 
geometries, with the position and size of extrema repro¬ 
duced well by both methods. For the two slit geometry, 
however, the conformal mapping method performs less 
well: the outer maxima are reproduced well, but the cen¬ 
tral one is not. This corresponds to the positioning of the 
plate where the edge of the pole piece coincides the cen¬ 
tral spar, where as discussed in section III B the mapping 



Figure 6: Comparison of finite element (solid lines) and con¬ 
formal mapping (points) approaches to calculating the local 
rate of dissipation P{q). 


is crowded. Clearly from this comparison, the method of 
choice appears to be finite elements, since the results are 
less noisy and all plate shapes can be simulated using the 
method. 


D. Effect on the Motion of the Pendulum 

Having formulated both models, we turn to the effect of 
the dissipation on the motion of the pendulum. For both 
the simplified wire model of section II and the magneto¬ 
quasistatic model of section III, the equation of motion 
has the form. 


X + (jJqX -I- \{x)x = 0. 
For the wire model. 




/ 3 - 


Bp 


L J p 


(39) 


(40) 


as read off from eq. (10); for the magneto-quasistatic 
model the Stokes drag is supplemented by the position- 
dependent dissipation rate, 

Xrngs = 2l3 + daPiq)Bl (41) 


As is well-known, the equation of motion (39) with 
constant A can be solved analytically yielding, 

x{t) = cos(w'f + 4>) + Xi (42) 


where xq and xi are amplitude and offset parameters, (j) 
is a phase and uj' = pijP — }?. Hence, the motion re¬ 
mains oscillatory for uj > X. If the plate is centered with 
respect to the pole pieces, xi = 0; qo and (j) similarly 
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I to power supply 


Figure 7: (a) Photograph of the experimental setup, designed 
as a classroom demonstration for introductory physics, (b) 
Detail showing placement of camera and screen for imaging, 
(c) Photograph of paddle shapes used. 

represent initial conditions that are experimentally con¬ 
trollable. Rather than fit to the entire trajectory, it is 
convenient to look only at the extreme points where the 
pendulum is stationary and cos (wot-!-(/') = ±1. Inserting 
this into (42) and rearranging, we obtain 

log \x{t)\= log xo - At/2, (43) 

showing that the values of the position of the pendulum 
at successive stationary points should obey a linear rela¬ 
tion on a plot of log(a;) versus t. 

For the magneto-quasistatic model, analytical solution 
of (39) is impossible and hence the equation must be 
solved numerically. This was performed in Mathematica. 
Because the strategy of looking at successive extrema on 
log scales remains a valuable fitting strategy even if the 
equation of motion is solved numerically, we addition¬ 
ally incorporated an extremum-finding routine to iden¬ 
tify the stationary points from the calculated trajectory. 
With these simplified analytical and numerical magneto¬ 
quasistatic models, we are now ready to compare them 
with experimental data. 


IV. EXPERIMENT 

In order to test the predictions of the calculations per¬ 
formed in previous sections, we performed an experiment 
to observe the motion of a pendulum with the different 
paddles attached and fit the models described in previous 
sections to the trajectories obtained. 

The magnetic brake apparatus used, shown in Figure 
7, was designed as a classroom demonstration for intro¬ 
ductory physics classes. The apparatus consists of an 
electromagnet with flat pole pieces of adjustable separa¬ 
tion, powered by a 12V car battery. A potential divider 
and ammeter were incorporated into the circuit so that 
the current through the magnet could be changed and 
measured. The pole pieces are 10 x 10cm and for the 


present experiment were separated by 2.5cm. The pen¬ 
dulum hangs from a supporting frame above the pole 
pieces and is made from a solid rod of length 39.5 cm 
to which different paddles of interest can be affixed with 
screws. Four plates were studied, each of side 11.43 cm 
by 11.43 cm with different arrangements of holes or slits 
as shown earlier in figure 1(d). 

Students collected videos of the motion of the pen¬ 
dulum using a cellphone camera, and in a later itera¬ 
tion a digital camera was used to achieve higher image 
quality and frame rate. Although data logging hardware 
could also be used where available, the video approach 
was adopted as it is inexpensive, required no additional 
apparatus beyond the existing demonstration and is eas¬ 
ily reproduced by others. A white piece of paper with 
a scale bar was attached to the pole piece behind the 
camera to provide a uniform background and facilitate 
quantitative analysis of the resulting videos. The cam¬ 
era was positioned so that the lens was aligned with the 
rod at rest; its position was marked with chalk so that it 
could be easily replaced after each trial. The plates were 
released from the same point outside of the magnet, de¬ 
termined by a ruler placed on top of the pole pieces. The 
plates were released when the camera started recording 
and data was collected until the pendulum came to a rest. 
Two sample frames from one of the movies are shown in 
figure 8(a). For each plate, 5 videos were collected with 
OA, 0.6A, lA, lAA and 1.8A current flowing through 
the circuit respectively. Selected videos were taken more 
than once to verify repeatability. These values obviously 
depend on the details of the circuit, but were chosen as 
follows: for the lower bound, a value was sought such that 
there was a visible difference between the motion of the 
plate with weakest damping (the one with 4 slits) in the 
on and off states. For the upper bound, we found a cur¬ 
rent level such that the plate with the strongest damping 
gave 2-3 oscillations before coming to a stop, i.e. we are 
in the oscillatory regime of w > A in equation (39) . 

Having collected the dataset of videos of the pendu¬ 
lum with different plates, the trajectories were extracted 
from the videos. This was done in two steps: First, the 
open-source image processing program Imaged [12] was 
used to “reslice” the movie, converting it from a stack of 
images of the (cc, y) plane at different time points t to 
a new stack of images of (x, t) at different vertical po¬ 
sitions y. An sample output is displayed in fig. 8(b), 
showing the damped motion of the pendulum as a func¬ 
tion of time for the row of the movie indicated by the 
dashed line in fig. 8(a). The second step in the data ex¬ 
traction was performed by a custom program written in 
Mathematica, included as Supplementary Material. This 
first subtracted the background, then for each time point 
located the position of the bob by finding the centroid of 
the pixels with intensity above a certain value. The result 
of this process was to yield for each movie a list of (x, t) 
representing the trajectory of the pendulum. Using the 
rulers and known frame rate of the camera (30fps), the 
list was calibrated to physical coordinates; the program 
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Figure 8: (a) Sample frames from a video of the pendulum’s 
motion. From this movie, the row of pixels indicated with 
the dashed line is taken from each frame and composited to 
produce the resliced image (b). 


also subtracts the equilibrium position of the pendulum 
from the x position. To facilitate comparison with the 
models, the extremum-finding routine mentioned in the 
previous section was used to identify a list of stationary 
points {(xiyti)} from each trajectory. A constant was 
subtracted from the times such that for each trajectory 
the first extremum was set to f = 0. 


V. MODEL COMPARISON 

Having obtained a set of experimental trajectories, we 
now fit both models to the data. Beginning with the 
simplified wire model, the analytical solution in section 
HID implies that for each individual trajectory the ex¬ 
trema should follow a straight line on semi-log scales with 
slope 

A„,.e = /3+f (44) 



r/r square 

/3 a 

Square 

1 

1.1 X 10"^ 0.060 

4 holes 

0.441 

4.4 X 10"® 0.130 

2 slits 

0.348 

1.2 X 10"^ 0.074 

4 slits 

0.136 

5.2 X 10"® 0.044 


Table I: Geometric parameters and fitted a, /3 for each plate 
using the simplified model. 


in principle be calculated or measured separately, but the 
purpose of fitting here is to establish whether the model is 
consistent with the results for all plates simultaneously. 
The slopes extracted from the ensemble of plates with 
different T = A^/L (calculated values are shown in table 
I) at different applied currents I should therefore univer¬ 
sally be described by, 

^wire — /3 oTJ^. (45) 

Estimates for the initial position of each trajectory xq 
were obtained by fitting a straight line to the extrema 
individually on semi-log scales. We then estimated /3 for 
each plate from the / = 0 trajectory for each plate; the 
results are shown in table I. Two pairs of plates seem 
to cluster around similar values; these were collected at 
different times indicating that the natural damping of 
the apparatus, perhaps lubricant in the hinge, changed 
between them. 

Results from the universal least-squares one-parameter 
fit for a are shown in fig. 9(a). As can be seen, the fit 
isn’t bad and follows many of the trajectories quite well. 
That said, there are many discrepancies. We therefore 
estimated a from each individual plate using the data at 
different applied field. While the details of these are not 
shown here, the quality of fit was very much better. The 
resultant values of a are displayed in table I. In particu¬ 
lar, it is not surprising that the plate with 4 holes is such 
an outlier as the wire model cannot account for this ge¬ 
ometry. Values for the other plates are the same order of 
magnitude, which together with fig. 9(a) demonstrates 
that the simplified model gives reasonable estimates for 
the influence of plate geometry even though it is based 
on quite egregious approximations. 

In the same spirit, we assume that the damping term 
in (39) for the magneto-quasistatic model can similarly 
be described by a universal constant a, 


From the set of results displayed as points in fig. 9, it 
is seen that this is generally the case, though some de¬ 
viations are apparent. In particular, the motion starts 
to decay more rapidly once \x/Lpp\ < 0.1; we therefore 
excluded this data from our fit. For each plate, several 
values of current / were used to take a succession of tra¬ 
jectories. Although it is obvious that Hq oc I, we did not 
explicitly measure the magnetic field and therefore ab¬ 
sorb this unknown constant of proportionality together 
with p the effective resistivity of the wire into a single fit¬ 
ting parameter a. This fitting parameter could, of course. 


>^mqs{x) =13 + aP{x)l'^. 

For this model, as well as xq, f3, a , P{x) and /, it is 
also necessary to specify the natural angular frequency 
w; the value from the 1 = 0 trajectory for each plate is 
used since [3 <+i ljJ for the weak damping in the absence 
of the magnetic field and hence uj' ~ ut. The results 
of the one-parameter universal fit for a are displayed in 
9(b). The quality of fit is significantly better, especially 
for large amplitudes and early times, and more consistent 
between plates than for the simple model, implying that 
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Figure 9: Fits for (a) the simplified wire model and (b) the magneto-quasistatic model. Experimental data is shown as points; 
the model is shown as lines. 


the magneto-quasistatic model yields good predictions of 
the effect of plate geometry. Of particular importance is 
the fact that the fit for the 4 hole plate is no worse than 
that of the other plates, suggesting that these features 
have been handled correctly by the model. 

Nonetheless, there remains an unexplained phe¬ 
nomenon in the data, i.e the anomalous increase in 
damping once the amplitude of the motion drops be¬ 
low x/Lpp ~ 0.1. This is not predicted by either model. 
The simplified model predicts a constant rate of damping, 
while the magneto-quasistatic model predicts a crossover 
from magnetically dominated damping at large ampli¬ 
tudes to Stokes dominated damping at small amplitudes. 
We tried adjusting the magneto-quasistatic model in var¬ 
ious ways—offsetting the position of the plate with re¬ 
spect to the pole pieces and changing the size of the plate 
relative to the pole pieces—but found no similar effect for 
any reasonable parameters for the experiment. Clearly, 
the Stokes form of the non-magnetic damping is most 
likely at fault, because the source of damping is likely to 
be primarily the hinge rather than air resistance; further 
analysis should be performed to determine the origin of 
this anomaly. 


VI. PEDAGOGICAL CONTEXT 

The scientific results presented in this paper were ob¬ 
tained in part by students in Spring 2013 as part of 
the Electromagnetic Theory II class at Tufts University, 
the second course in the E&M sequence in the gradu¬ 
ate program at Tufts; two students from the class are 
co-authors on the paper (CW and JW) having spent a 
considerable time subsequently acquiring additional data 
and performing calculations. In this section, the peda¬ 
gogical context of these activities is briefly documented 
in the hope that it is of use to others. The idea of incor¬ 


porating (simple) experiments in class was motivated by 
an observation that key scientific skills expected of gradu¬ 
ate students—connecting theory and experiment; article 
reading; journal selection and article preparation—are 
not included explicitly in the curriculum. The project 
was spread over the second half of the semester as fol¬ 
lows:— 

1. A prepatory homework where students derived the 
magneto-quasistatic model III. 

2. Two Ij hour class periods where students worked 
in groups. The students were divided into groups 
of four and worked on: i) conformal mapping anal¬ 
ysis, ii) finite element analysis and iii) experiments. 
The instructor provided scaffolding activities to 
help students learn FlexPDE and the Schwarz- 
Christoffel mapping technique. These activities 
continued subsequently outside of class. 

3. A reading activity whereby students had to identify 
possible journals for publication of this work; this 
was followed by a 15 minute discussion of journal 
selection in class. 

4. A Just-In-Time Teaching[13]-like pre-class activity 
in which students had to identify the structural el¬ 
ements in a typical paper and construct an outline; 
these were passed round anonymously for peer feed¬ 
back. 

5. A homework activity in which students prepared 
an introduction for the paper. Students were asked 
to read each other’s work and identify strong fea¬ 
tures as well as possible modifications as a pre-class 
activity. This was followed by a discussion in class. 

6. Each group collaboratively produced a draft sec¬ 
tion for the final paper as part of a homework; these 
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were combined by TJA into a coherent document 
and extensively rewritten during the second itera¬ 
tion of experiments. 

The focus of the class and restricted time available neces¬ 
sitated some design trade-offs. For example, we chose to 
use an existing commercial software package, FlexPDE, 
as opposed to implementing a custom finite element pro¬ 
gram, for several reasons: First, the scale and com¬ 
plexity of problems solvable with a specialist code is 
much greater than with a simple custom-written pro¬ 
gram and hence should be applicable to situations en¬ 
countered in student’s research projects. Second, the 
program contains advanced features such as adaptive re¬ 
finement—allocation of grid points based on estimated 
error—that are time-consuming to implement. Third, by 
removing the focus of the exercise from programmatic im¬ 
plementation, time can be spent on understanding how to 
properly interpret the output using fundamental princi¬ 
ples of numerical analysis such as order, error estimation 
and convergence. 

While the project-based approach is well-grounded in 
a theory of learning (constructivism), due to the small 
= 13 size of the class, it is difficult to rigorously 
demonstrate the effectiveness of the activity. Future it¬ 
erations would strongly benefit from pre and post test¬ 
ing. Even so, the outcomes of the project indicate that 
students made considerable progress towards learning 
these challenging scientific skills. Moreover, the class was 
very highly rated in the formal feedback mechanism. It 
is hoped that, in future, physics education researchers 
might attempt a thorough analysis of this and related 
approaches at the graduate level. 


VII. CONCLUSION 

The present work has provided a careful, quantitative 
analysis of a familiar classroom demonstration, the van 
Waltenhofen pendulum. A magneto-quasistatic model 
has been shown to successfully predict the relative damp¬ 
ing of a selection of plates with different geometry using 
conformal mapping and finite elements as tools to carry 
out the necessary calculations. As a by-product, a num¬ 
ber of valuable resources for more general use, including 
plots of the current distribution during the course of the 
motion of the pendulum, have been produced. Addition¬ 
ally, a greatly simplified model, suitable for use in intro¬ 
ductory physics classes, has been formulated that could 
be used to make the classroom demonstration more quan¬ 
titative. 
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